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ABSTRACT: We explore different branches of the fermion doublers with Wilson fermion in per- 
turbation theory, in the context of additive mass renormalization and chiral anomaly, and show that 
by appropriately averaging over suitably chosen branches one can reduce cut-off artifacts. Com- 
paring the central branch with all other branches, we find that the central branch, among all the 
avatars of the Wilson fermion, is the most suitable candidate for exploring near conformal lattice 
field theories. 



Contents 



1 Introduction 



1 



2 Preliminaries 



2 



3 Additive Renormalization in Fermion Self Energy 



2 



4 Chiral Anomaly 



4 



5 Discussion and Conclusions 



6 



1 Introduction 

The most reliable techniques to investigate various non-perturbative aspects of quantum field the- 
ories are provided by lattice methods. Putting fermion on a lattice, however, has turned out to be 
highly non-trivial because of the notorious doubling problem. Naive discretization of the Dirac 
action leads to 16 solutions (called doublers) in the four dimensional theory. Among the various 
solutions suggested to cure this problem, Wilson fermions [1] are conceptually the simplest and 
straightforward to implement. In the Wilson approach a dimension five operator is added to the 
action thereby sending the masses of the extra fifteen fermions to the order of the cutoff. Thus the 
extra fermions decouple in the continuum limit. It is well-known that the Wilson formulation of 
lattice gauge theory preserves discrete symmetries of the continuum formulation which simplifies 
the construction of lattice operators that correspond to the observables in the continuum theory. 
The Wilson term in the action, which is introduced to remove the doublers, breaks chiral symmetry 
in accordance with the well-known Nielson-Ninomiya no-go theorem [2]. Wilson term reproduces 
the correct axial anomaly [3-5] eventhough it leads to additive renormalization for the fermion 
mass. 

The sixteen doublers are classified into five branches. Almost all of the studies so far, both 
analytical and numerical, have focused on the so-called first (physical) branch. However, very 
recently, occurrence of an enhanced symmetry has been discovered in the central branch [6-8] 
when the on-site terms (mass term and that from the Wilson term) are absent in the action. The 
enhanced symmetry prohibits additive renormalization through radiative corrections. Since in this 
case, the central branch yields six massless fermions, as suggested by ref. [7], an alternative way 
to simulate twelve flavour non-abelian gauge theories emerges. Such theories are of interest in the 
context of beyond standard model physics (for recent reviews, see [9-12]). 

In this work, by introducing a generalized Wilson term containing a branch selector index, 
we investigate the additive fermion mass shift and chiral anomaly to 0{g ) in lattice perturbation 
theory for all the branches of the fermion doublers. 
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2 Preliminaries 



We denote the generalized Wilson fermion action by 

S F [\ir,\jr,U](i B ) = a A £^M*y(ia)yr y = a 4 J^W* [j l iD ll + W{i B )+m] w y with (2.1) 



£^r,/i ^e+jUj U x _pn 8x-H,y 



and 



r T 1 



(2.2) 
(2.3) 



The factor ;'b is the branch selector index which takes the values 0, 1,2, 3, and 4 for first, second, 
third (central), fourth and fifth branch of the doubler respectively. We consider the transformations 
[7] 



(2.4) 



The action is invariant under these transformations but for the local terms. For m = and = 2 
(massless limit of the central branch), the action thus possesses this additional symmetry which 
prevents additive renormalization of the fermion mass through radiative corrections. 



3 Additive Renormalization in Fermion Self Energy 

In this section we calculate the additive shift to G (g 2 ) in the fermion mass (for am = 0) using 
lattice perturbation theoiy [13]. 

The additive shift in fermion mass due to the tadpole diagram (figure 1) is 

r 1 

8m = - - -g 2 C F Ycos(apn) Z (3.1) 
a 2 ~ 

/ \ _1 

d 4 k | at „:„2fak k \\ „„j ^ N 2 -\ 



with Z Q = J^\4U sin 2 ( ^ ) ] and C F = ^ for SU(N) 




Figure 1. Tadpole diagram 
Results for different branches are as follows. 
First branch: ap^ = (0, 0, 0, 0) — > 8m = —2^g 2 Cp Zq . 

Second branch: ap^ = (n, 0, 0, 0) or any of the other three permutations — > 8m = — l^g 2 Cp Zq . 
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Third (central) branch: ap^ = (ft, ft, 0, 0) or any of the other five permutations — > 8m = 0. 
Fourth branch: ap^ = (ft, ft, ft, 0) or any of the other three permutations — > 8m = +l^g 2 Cp Z ( 



o • 



Fifth branch: ap^ = (ft, ft, 71, ft) — > 8m = +2-g Cf Z ( 
Next consider the additive mass shift in fermion mass due to sunset diagram. The gauge boson 
propagator in Feynman gauge is given by 

^^^{^Isin 2 ^^}" 1 = S MV S«* {(1/a 2 ) W p ^y l , (3.2) 
whereas the fermion propagator has the form 



with 



S lm (k,i B ) = 8"4 £ iy „™M + -J( k {i B ) 
I v- a a 



and the fermion-gauge boson vertex is 

tnr(M = -s(n { irp co S ?^ + ™„?^} . 

Then the fermion self energy from sunset diagram can be evaluated as 



d 4 k 



(3.3) 



(3.4) 



(3.5) 



(3.6) 



The additive mass shift arising from the fermion self energy (sunset) can be written as 

d 4 k N r 



8m = -g Cp 
a 



(2ft) 4 D r 



(3.7) 



where D, = W p>k (F 2 + r 2 ^C 2 (i B )^ with F 2 = £^ sin 2 (ak^). We introduce = sin(a^), 
S p = sin(^£-) and C p = cos(^-). The expressions for N r and D r for different branches are given 
below. 

First branch: ap^ = (0, 0, 0, 0). 



4 r 



N r = £ J( k (is = 0) (r 2 S 2 - C 2 ) + r 
P =i 



D r = Wr 



p.k 



F 1 + r 1 ^{i B = Q) 



(3.8) 
(3.9) 



Second branch: ap^ = (ft, 0, 0, 0) or three other permutations. Explicitly for ap^ = (ft, 0, 0, 0) 



Nr 



Jt k {i B = l){S 2 p -r 2 C 2 p )+F 2 p ] _ +£ \^ k (i B = l)(r 2 S 2 p -C 2 p )+F i 



p=l 



p=2 



D r = W„ k 



F 1 + r 2 Jtj L {i B = \) 



(3.10) 
(3.11) 
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Figure 2. Sunset diagram 



Third (central) branch: ap^ = (71, 71, 0, 0) or five other permutations. 
Explicitly for ap^ = (n, n, 0, 0) 

2 4 

N r = - £ Jt k (i B = 2) {S 2 p - r 2 C 2 p ) + r 2 p + £ J? k {i B = 2) (r 2 S 2 p - C 2 , ) + P 

p=l P=3 L 



p.k 



l (iB = 2) 



(3-12) 



(3.13) 



Fourth branch: ap^ = (n, n, n, 0) or three other permutations. Explicitly for ap^ = (n, 7T, 7C, 0) 



N r 



I 



Ui B = 3)(S 2 p -r 2 C 2 p ) + r 2 p + ^ k (i B = 3)(r 2 S 2 p -C 2 p ) + r] 



.,2 S 2 



p=4 



, (3.14) 



Fifth branch: ap^ 



r 2 + r 2 ^f (i fl = 3) 
= (7T, 7T, 7T, 7T). 



AT, 



4 

I 

p=l 



«* = 4)(S 



- 2 ^)+r p 



D r = Wr 



p.k 



T 2 + r 2 Jtj;(i B = 4) 



(3.15) 



(3.16) 



(3.17) 



In table 1 we present the numerical values of the additive mass shift separately from sunset 
and tadpole contributions for the fermion at different branches. In figure 3 we plot the magnitude 
of the total additive mass shift (tadpole + sunset) versus the branch number. Note that the shift is 
maximum for the first and the fifth branches and is minimum (zero) for the third (central) branch. 
The absence of additive renormalization in fermion self energy for the central branch to &{g 2 ) is 
explicitly shown also in ref. [8]. 



4 Chiral Anomaly 

Now we consider the flavor singlet axial Ward Identity 

(Aps^x)) = 2m{\ff x Y 5 Yx) + (Xx) 



(4.1) 
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Drdncn 


8 m/(g 2 C F ) 


Sunset 


Tadpole 


First 


-0.0158 


-0.3099 


Second 


+0.0148 


-0.1549 


Third 


0.0000 


0.0000 


Fourth 


-0.0148 


+0.1549 


Fifth 


+0.0158 


+0.3099 



Table 1. Numerical values of the additive mass shift for fermion at different branches for r = 1 and L = 200. 
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Figure 3. The magnitude of the total additive mass shift (tadpole + sunset) plotted versus the branch number. 



where (&) denotes the functional average of G. Explanation of other terms are as follows: 
1 



Aj/(x) 
and (&) 



Here 



[f(x) — f(x — ju)] , the backward derivative, 



WxYn 7s Ux lt i Wx+pl + Wx+ii 7m 7s Wx 

1 



the axial vector cuiTent 



(4.2) 
(4.3) 



2 g 2 e^px Trace F^ v (x) F pX {x) ——j £ cos(p^a) cos(p v a) cos(p p a) 



x r 



„2 



(2^) 4 t 

cos(p A a) [am + r^r p (i B )]- 4rsm 2 (p x a) {%{i B )f , 



£ ^vpA Trace fjivW f pA W I(am,r,L). 



(4.4) 



f^sin 2 (ap M ) + 


ara + r. 


£ p {i B ) 













(4.5) 



-5- 



Explicitly, £ p = (^) 4 Eni,« 2 «3 «4 - m an our plots it is the anomaly integral denoted by the function 
I(am,r,L) which we have plotted. 

Following Karsten and Smit [3], the limits on the momentum sum are changed from (— n, +7l) 
to (— 7r/2,37r/2) and further the momentum sum hypercube is divided into 16 smaller hypercubes 
corresponding to (— lt/2, +7l/2) and (+7T/2, +37T/2) for each ap^,/j. = 1,2,3,4. Thus the total 
anomaly contribution is decomposed into the contributions from five different types of species and 
the anomaly integral takes the form / = To — 4I\ + 6h — ^h + h- In all the four momentum 
integrations range from (—n/2, +7T/2) and in h they range from (+7T/2, +37T/2). In I\ one of the 
momentum integrations range from (+7T/2, +37T/2), the rest from (— 7r/2, +7i/2) and vice-versa 
for I3. In the third (central) branch h two momentum integrations range from (+7T/2, +37T/2) and 
the rest from (— 7T/2, +7i/2). 

First, to perform the integration analytically, we set the bare mass am = 0, use the identity [3] 



^pih) cos(app) — 4r^# p (/g) sin 2 (app) 



d(ap p ) 



1 v 1 
^pih) +2 - sin ( a Po) 



sinfap^K J^p{iB) +^sin 2 (£?p f7 



(4.6) 



and do a partial integration. In the infinite volume continuum limit, the results for the integrals are 
as follows. 

First branch: Iq — > 1, I\, I2, h, h — > 0. 
Second branch: 1\ — > 1, Iq, h, h, h — > 0. 
Third (central) branch: I2 — > 1, Io, hi hi h -^-0. 
Fourth branch: h — > 1> lo, h, h, h —^0. 
Fifth branch: I4 — >• 1, 7 , h, h, h — > 0. 

Since numerical simulations are performed at finite volume and finite lattice spacing, it is of interest 
to study the effect of symmetry violation on the anomaly integral as a function of the lattice fermion 
mass at finite volume [14, 15]. In order to avoid the zero mode problem we have used antiperiodic 
boundary condition in all four directions. In figure 4 we plot the function I(am,r,L) for r = 1.0 
and L = 100 as a function of am for the first and fifth branches (left) and for the second and fourth 
branches (right). In figure 5, we plot the function I(am,r,L) for r = 1.0 and L = 100 as a function 
of am for the central branch. From figure 4 (left), we observe that the cut-off effects are almost 
equal and opposite for first and fifth branches. Similar observation can be made regarding second 
and fourth branches from figure 4 (right). Comparing figures 4 and 5, we conclude that cut-off 
effects are minimal for the central branch. We have picked L = 100 for our plots as we have found 
that finite volume effects are negligible at this volume for the range of am shown in the figures. 



5 Discussion and Conclusions 

It is well known that the naive discretization of the fermionic action gives rise to sixteen degenerate 
species including the desired physical one. These sixteen species are grouped into five branches 
with degeneracy (chirality) given by 1(1), 4(-l), 6(1), 4(-l) and 1(1), rendering the theory free of 
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Figure 4. The function I(am, r,L) for r = 1 .0 and L = 100 as a function of am for the first and fifth branches 
(left) and for the second and fourth branches (right). 
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Figure 5. The function I(am, r,L) for the central branch for r = 1 .0 and L = 100 as a function of am. 



chiral anomaly. With the conventional Wilson term in the continuum limit, apart from the first 
branch, species corresponding to all other branches become infinitely massive and decouple from 
the theory thereby reproducing the correct chiral anomaly. The branches other than the first one are 
rarely explored. However, recently the existence of an additional symmetry in the central branch 
which prohibits additive renormalization of fermion mass has been discovered in the ref. [7]. 

In this work, in order to explore all branches we introduce a generalized Wilson term con- 
taining a branch selector index By choosing one can make the fermions belonging to a 
particular branch physical. The fermions belonging to the rest of the branches become infinitely 
massive and decouple from the theoiy in the continuum limit. The conventional Wilson term cor- 
responds to is = 0. To investigate the effect of radiative corrections, we calculate the additive mass 



renormalization in fermion self-energy and the chiral anomaly to 0(g 2 ) in perturbation theory for 
all the branches. 

First we summarize the results of additive mass shift from tadpole and sunset contributions. 
The tadpole contributions for the first and fifth branches are equal in magnitude but opposite in 
sign. Same is true for the sunset contributions also. Thus 8m vanishes if we average over the first 
and fifth branches. Similar observations hold for the second and the fourth branches also. Coming 
to the central branch the additive mass shifts from tadpole and sunset contributions separately 
vanish. This leads to the absence of additive mass renormalization in accordance with theoretical 
expectation. In the calculation of chiral anomaly first we perform an analytical calculation setting 
am = and using the Karsten-Smit identity. We find the correct value of the anomaly for different 
branches with corresponding degeneracy factors and signs dictated by the chiral charges in the 
continuum limit. Since numerical simulations are performed at finite volume, finite lattice spacing 
and finite fermion mass, we have studied the effect of symmetry violation (given in eq. 2.4) on the 
anomaly integral as a function of the lattice fermion mass. The cut-off effects are almost equal in 
magnitude but opposite in sign for the first and the fifth branches. Same holds for the second and 
the fourth branches also. The cut-off effect is minimal for the central branch. 

In conclusion, our exploration of the different branches of the fermion doublers in perturbation 
theory, in the context of additive mass renormalization and chiral anomaly, has shown that by 
appropriately averaging over suitably chosen branches one can reduce cut-off artifacts. Comparing 
the central branch with all other branches, we find that the central branch, among all the avatars 
of the Wilson fermion, is the most suitable candidate for exploring near conformal lattice field 
theories. 
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